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Abstract 

Recent work in low energy pion physics is reviewed. One of the exciting 
new developments in this field is that simulations of QCD on a lattice now 
start providing information about the low energy structure of the contin- 
uum theory, for physical values of the quark masses. Although the various 
sources of systematic error yet need to be explored more thoroughly, the 
results obtained for the correlation function of the axial current with the 
quantum numbers of the pion already have important implications for the 
effective Lagrangian of QCD. The consequences for mr scattering are dis- 
cussed in some detail. The second part of the report briefly reviews recent 
developments in the dispersion theory of the scattering amplitude. One of 
the important results here is that the position of the lowest resonances of 
QCD can now be determined in a model independent manner and rather 
precisely. Beyond any doubt, the partial wave with I = £ = contains a 
pole on the second sheet, not far from the threshold: the lowest resonance 
of QCD carries the quantum numbers of the vacuum. 
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1 Introduction 

The pions are the lightest hadrons and we know why they are so light. Since 
the underlying approximate symmetry also determines their basic properties at 
low energy, the interaction among the pions is understood very well. In fact, in 
the threshold region, the 7T7t scattering amplitude is now known to an amazing 
degree of accuracy. pQ In particular, we know how to calculate the mass and 
width of the lowest resonance of QCD in a model independent manner. [3j The 
actual uncertainty in the pole position is smaller than the estimate given in 
the 2006 edition of the Review of Particle Physics, [3] by more than an order of 
magnitude. 

The progress made in this field heavily relies on the fact that the dispersion 
theory of 7T7t scattering is particularly simple: the s-, t- and u-channels represent 
the same physical process. As a consequence, the scattering amplitude can be 
represented as a dispersion integral over the imaginary part and the integral 
exclusively extends over the physical region. [1 Throughout the following, I work 
in the isospin limit (m u = to<j = m and e = 0), where the representation 



involves only two subtraction constants These may be identified with the S- 
wave scattering lengths <io> a o- The projection of the amplitude on the partial 
waves leads to a dispersive representation for these, the Roy equations. For a 
thorough discussion, I refer to ACGL.JS] 

The pioneering work on the physics of the Roy equations was carried out 
more than 30 years ago. [5] The main problem encountered at that time was that 
the two subtraction constants occurring in these equations were not known: if 
the values of (Zq, Oq are contained in the so-called universal band - the region 
spanned by the two thick lines in figure 1 below - the Roy equations admit 
a solution. Since the data available at the time were consistent with a very 
broad range of 5-wave scattering lengths, the Roy equation analysis was not 
conclusive. 

2 Low energy theorems for the /S-wave scatter- 
ing lengths 

The insights gained by means of chiral perturbation theory (%PT ) thoroughly 
changed the situation. [7] The corrections to Weinberg's low energy theorems 8 
for a[], <Zq (left dot in figure 1) have been worked out to first non- leading order [9] 
(middle dot) and those of next-to-next-to leading order are also known[10] (dot 
on the right). As demonstrated in CGL,[1] the chiral perturbation series con- 
verges particularly rapidly near the center of the Mandelstam triangle, so that 
very accurate predictions for the scattering lengths are obtained by matching 
the chiral and dispersive representations there. Using this method, the low 
energy theorems for <Zq, a§ may be brought to the form 

2a» + 7a= = - 6t (t, - gjj) + Mja, + O(Mj) , (la) 

2aS _ 5o ; . g| + 24 ,(^) 4 (^gZ) + M. a4 + 0(M;) ^ lt) 

While the chiral expansion of a^a^ starts at order M%, with a term that is 
fixed by the pion decay constant, the combination 2a[j + 7oq vanishes at leading 
order. Two types of corrections occur at O(M^): a contribution from the cou- 
pling constants £3, £4 of the effective Lagrangian and one involving dispersion 
integrals, which I denote by M^a^ and M^a^, respectiveljH (these also contain 
higher order contributions). 

A representation similar to (1) was given earlier, [9] but there, the disper- 
sion integrals M^a-i, M^oti were evaluated at leading order, where they can be 
expressed in terms of the D-wave scattering lengths. The virtue of the above 

The value used for the basic QCD parameters in this theoretical limit is a matter of 
convention. It is convenient to choose K.qqd an d m such that i*V and M-k agree with the 
observed values of F + and M + . The result quoted by the PDG|3] corresponds to F + = 
92.42 ± 0.10 ± 0.26 MeV. The scattering lengths are given in units of M Jr+ . 
2 In the notation of CGL: 03 = 2ao + 7«2, 04 = 2ao — 5a2- 



2 




Figure 1: S wave scattering lengths: theoretical results. 

form of the low energy theorems is that (a) the neglected terms of order M® are 
much smaller than in the straightforward chiral perturbation series and (b) the 
dispersion integrals can be evaluated quite accurately, on the basis of the Roy 
equations. Since the integrals converge very rapidly, they depend almost exclu- 
sively on the scattering lengths a[J, ag. In the vicinity of the physical values, the 
linear representations 

M*a 3 = 0.135 + 0.77 (ag - 0.220)- 1.50 (a 2 + 0.0444), (2a) 
Af*a 4 = 0.061 + 0.48 (a° - 0.220) - 0.26 (ag + 0.0444) , (26) 

provide an excellent approximation. For a detailed discussion, in particular also 
of the uncertainties, I refer to CGL. The estimates given there show that the 
contributions of order M® in (la) generate a correction of order 0.002, while 
those in equation (lb) are of order 0.003. 

The low energy theorems (1) show that the S- wave scattering lengths are 
related to the coupling constants £3 and £4. Indeed, these play a central role 
in the effective theory, because they determine the dependence of M w and F v 
on the quark mass at first non-leading order: the expansion of M 2 and F v in 
powers of the quark mass starts with 

Ml = AI 2 {1 + ±x£ 3 + 0(x 2 )} , F^ = F{l + xh + 0(x 2 )} , (3) 

where M 2 = 2Bm is proportional to the quark mass, F is the value of the pion 
decay constant in the chiral limit, and x = M 2 / (16ir 2 F 2 ). A crude estimate 
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for £3 can be obtained from the mass spectrum of the pseudoscalar octet: £3 = 
2.9 ± 2.4. [9] For £4, an analogous estimate (£4 — 4.3 ± 0.9) follows from the 
experimental value of the ratio Fk/F^. The low energy theorem for the radius 
of the scalar pion form factor yields a more accurate result: £4 = 4.4 ± 0.2. [TT] 
The error is smaller here, because the theorem holds within SU(2)xSU(2), so 
that an expansion in the mass of the strange quark is not necessary. The small 
ellipse in figure [T] shows the prediction for the scattering lengths obtained in 
CGL with these values of £3 and £4 (since the residual errors are small, the 
corrections of order M® are not entirely negligible - the curve shown includes 
an estimate for these). 

The narrow strip that runs near the lower edge of the universal band indicates 
the region allowed if the coupling constant £3 is treated as a free parameter, while 
£4 is fixed with the scalar radius, like for the ellipse. By construction, the upper 
and lower edges of the strip are tangent to the ellipse. As implied by equations 
(lb) and (2b), the strip imposes an approximately linear correlation between 
the two scattering lengths - a curvature becomes visible only in the region where 
the theoretical prediction would be totally wrong. 

3 Lattice results 

On the lattice, dynamical quarks can now be made sufficiently light to establish 
contact with the effective low energy theory of QCD. In particular, the MILC 
collaboration obtained an estimate for the coupling constants L4, L5, L§. Lg, 
of the effective chiral SU(3)xSU(3) Lagrangian. [T2"] Using standard one loop 
formulae. [T3] these results lead to £3 = 0.8 ±2.3 and £4 — 4.0 ±0.6. In the plane 
of the two S-w&ve scattering lengths, the MILC results select the region indi- 
cated by the second ellipse shown in figure [TJ which is slightly larger. It would 
be worthwhile to analyze these data within SU(2)xSU(2), in order to extract 
the constants £3 and £4 more directly. The one loop formulae of SU(3)xSU(3) 
are subject to inherently larger corrections from higher orders, so that the re- 
sults necessarily come with a larger error (note that, at the present state of the 
art, the analysis of the lattice data relies on the one loop formulae). 

The second narrow strip shown in the figure is obtained by treating £4 as 
an unknown. Instead of showing the region that belongs to the value of £3 
used for the ellipse, I am taking £3 from a lattice study, which appeared very 
recently. [2] In that work, the dependence of M% and F^ on the quark mass 
is investigated for two flavours of light Wilson quarks. The results for M% are 
consistent with one loop xPT . With M Klaf = 495 MeV, the result £3 = 0.5(5) (1) 
implies £3 — 3.0 ± 0.5 ± 0.1. Since this is close to the center of the estimated 
range, the strip runs through the middle of the ellipse. Taken at face value, 
the uncertainty in the lattice result for £3 is four times smaller than the one 
in the estimate £3 — 2.9 ± 2.4, obtained from the SU(3) mass formulae for 
the pseudoscalars, more than 20 years ago. [9] While the width of the ellipse is 
dominated by the uncertainty in the value used for £3, those associated with 
the higher order contributions and with the phenomenological uncertainties in 
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the dispersion integrals do affect the width of the strip belonging to the value 
of Del Debbio et al. - this is why the higher precision for £3 reduces the width 
by less than a factor of four. 

The coupling constants depend logarithmically on the pion mass: 

so that their value changes if the quark mass is varied. The quoted value for £3 
corresponds to A 3 ~ 560 MeV. For < A3, the curvature term is positive and 
reaches a maximum at M v = A^/y/e ~ 340 MeV. Even there, \xt$ is less than 
0.05, so that M%/m is indeed nearly a constant. Although more work is needed 
to clarify all sources of uncertainty, the calculation beautifully confirms that, to 
a very good approximation, M% is linear in the quark mass. Hence the quark 
condensate indeed represents the leading order parameter of the spontaneously 
broken symmetry. |15j The same lattice data also yield information about the 
dependence of on the quark mass, but the comparison with the one loop 
formula of xPT is not conclusive in this case: the result for £4 depends on how 
the data are analyzed. [14] 

The quark mass dependence of M w and is also being studied by the 
European Twisted Mass Collaboration. The preliminary results for the relevant 
effective coupling constants read £3 = 3.5 and £4 — 4.4. A thorough analysis of 
the uncertainties is under way. [IB] 

Finally, I mention that the scattering length of the exotic S-wave, <Zq, can 
be determined directly, from the volume dependence of the energy levels occur- 
ring on the lattice. The horizontal band shown in figure Q] indicates the result 
obtained in this way by NPLQCD.p~7] It is also consistent with the other pieces 
of information shown in the figure. Although possible in principle, it is difficult 
to extend this method to the isoscalar channel. 



4 Experiment 

Since the pions are not stable, they first need to be produced before they can be 
studied, so that the experimental information is of limited accuracy. Moreover, 
there are inconsistencies among the various data sets. One of the problems 
is that, in most production experiments, the two pions in the final state are 
accompanied by other hadrons. I know of three exceptions: production via 
photons in e + e~ collisions or via VF-bosons in the decays — > vit^it® and 
K — > eviTTT. The first two yield excellent information about the electromagnetic 
and weak vector form factors of the pion and hence also about the P-wave 
7T7T phase shift. The electromagnetic form factor is of particular interest in 
view of the Standard Model prediction for the magnetic moment of the muon. 
In this connection, the theoretical understanding of the final state interaction 
among the pions achieved in recent years can be used to reduce the experimental 
uncertainties, in particular in the region below 600 MeV, where the data are 
meager. [18] The K to tttt transition form factors relevant for the third process 
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Figure 2: S wave scattering lengths. Experimental results 

allow a measurement of the phase difference between the S- and P-waves, <5g— S\. 
The ellipse labeled E865[19. in figure 2 shows the constraint imposed on a® and 
Oq by the K e 4 data collected at Brookhaven. The two ellipses denoted DFGS[20] 
represent the la and 2a contours obtained by combining these with other irir 
data. The region labeled PY[5T] has roughly the same experimental basis. 

A very interesting proposal due to Cabibbo|22] has recently been explored 
in the NA48/2 experiment at CERN: the cusp seen near threshold in the decay 
K ± — ► ^ir ^ allows a measurement of the difference <Zq — Oq of scattering 
lengths. [53] In this case, the final state involves three pions, so that the analysis 
of the data is not a trivial matter. pMl [231 [23] As seen in the figure, the result 
is in good agreement with the Brookhaven data, as well as with the theoretical 
predictions. The same group has also investigated a very large sample of K e 4 
decays. [27] The result for the scattering length Og is indicated by the broad 
vertical band. It is not in good agreement with E865, nor with their K — > 3n 
result. The discrepancy calls for clarification. 

An ideal laboratory for exploring the low energy properties of the pions is 
the atom consisting of a pair of charged pions, also referred to as pionium. 
The DIRAC collaboration at CERN has demonstrated that it is possible to 
generate such atoms and to measure their lifetime. [55] Since the physics of the 
bound state is well understood. [23] there is a sharp theoretical prediction for 
the lifetime. [T] The band labeled DIRAC in the figure shows that the observed 
lifetime confirms this prediction. Pionium level splittings would offer a clean 
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and direct measurement of the second subtraction constant. Data on irK atoms 
would also be very valuable, as they would allow to explore the role played by 
the strange quarks in the QCD vacuum. [551 130] 



5 Dispersion relations 

After the above extensive discussion of our current knowledge of the two subtrac- 
tion constants, I now wish to discuss their role in the low energy analysis of the 
7T7T scattering amplitude, avoiding technical machinery as much as possible. Al- 
though the Roy equations represent an optimal and comprehensive framework 
for that analysis, the main points can be seen in a simpler context: forward 
dispersion relations. [31] More specifically, I consider the component of the scat- 
tering amplitude with s-channel isospin I = 0, which I denote by T°(s,t). It 
satisfies a twice subtracted fixed- i dispersion relation in the variable s. In the 
forward direction, t — 0, this relation reads 

n, \ s(s -AMI) f°° dxImT°(x,0) 

Rel>,0) . C0 + Cl , + J—^ P J^ x{x _ iM ^^ s) + (5) 

s{s-AM%) f°° dx{lmT°(x,0)-3ImT 1 (x,0)+5lmT 2 (x,0)} 
+ n 7 4M 2 3x{x-AM%) (x + s - AM%) ' 

The symbol P indicates that the principal value must be taken. The first inte- 
gral accounts for the discontinuity across the right hand cut, while the second 
represents the analogous contribution from the left hand cut, where the com- 
ponents of the scattering amplitude with I = 1, 2 also show up. According to 
the optical theorem, the imaginary part of the forward scattering amplitude 
represents the total cross section: in the normalization of ACGL,[5] we have 
ImT 7 (s,0) = \/s(s — AM%) a( ot (s). In this notation, the physical total cross 
sections are given by 

a tot - G tof a tot - 2 a tot+ 2 a tot> \ b ) 

Tr^TT* _ 10 i 1-1 I 1-2 rr^"^" — I/T° _]_ Irr 2 

"tot — 3 "tot "i 2 tot "r" 6 tat ' "tot ~ 3° tot ' 3" tot ■ 

As mentioned already, the subtraction term is determined by a®, a 2 ,: 

co +a s = 32n |a° + (2a° 5a 2 ) -J^qfj ■ (7) 

A dispersion relation of the above type also holds for other processes. What 
is special about rnr is that the contribution from the crossed channels can be 
expressed in terms of observable quantities - total cross sections in the case of 
forward scattering. The contribution from the left hand cut is dominated by 
the p-meson, which generates a pronounced peak in the total cross section with 
1 = 1. This contribution is known very accurately from the process e + e~ — > 
7r + 7r~ . In the physical region, s > AM%, the entire contribution from the crossed 
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channels is a smooth function that varies only slowly with the energy. Note, 
however, that this contribution is by no means small. [32] 

The angular momentum barrier suppresses the higher partial waves: at low 
energies, the first term in the partial wave decomposition 

ReT°(s,t) = 327r{Reig(s)+5P 2 (z)Re^(s) + ...} (8) 
* = $(4M*-8)(l-z) 

represents the most important contribution. In the vicinity of the threshold, 
where the contribution from the higher angular momenta is negligibly small, 
the dispersion relation ((5]) thus amounts to an expression for the real part of 
the isoscalar S- wave. For brevity, I refer to this partial wave as S°. 

Equation ([5]) imposes a very strong constraint on S°, for the following reason. 
Suppose that all partial waves except this one are known. The relation then 
determines Retg(s) as an integral over Im^(i). In the elastic region, unitarity 
already fixes the real part in terms of the imaginary part - we thus have two 
equations for the two unknowns Re£g(s) an d Imi[j(s). Accordingly, we may 
expect that, in the elastic region, the f.d.r. unambiguously fixes the S^-wave. 
Indeed, this is borne out by the calculation, which is described in some detail 
elsewhere. [33] 

One of the main results established in CGL is that the Roy equations fix 
the behaviour of the S°-wave below 800 MeV almost entirely in terms of three 
parameters: the two subtraction constants ao,flo an( ^ the- value of the phase 
at 800 MeV. In particular, the behaviour of the scattering amplitude at high 
energies is not important. This can also be seen on the basis of equation ([5]): 
replacing our representation of the scattering amplitude above KK threshold 
by the one proposed in KPY.[31] for instance, the solution of the f.d.r. very 
closely follows the solution of the Roy equations that belongs to the same value 
of the phase at 800 MeV. [33] 

6 The lowest resonances of QCD 

The positions of the poles in the S-matrix represent universal properties of the 
strong interaction, which are unambiguous even if the width of the correspond- 
ing resonance turns out to be large, [M] but they concern the non-perturbative 
domain, where an analysis in terms of the local degrees of freedom of QCD - 
quarks and gluons - is not in sight. First quenched lattice explorations of the 
pole from the a have appeared. [35] but in view of the strong final state inter- 
action in this channel, it will take some time before this state can reliably be 
reached on the lattice. 

One of the reasons why the values for the pole position of the a quoted by 
the Particle Data Group cover a very broad range is that all of these rely on the 
extrapolation of hand made parametrizations: the data are represented in terms 
of suitable functions on the real axis and the position of the pole is determined 
by continuing this representation into the complex plane. If the width of the 
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resonance is small, the ambiguities inherent in the choice of the paramctrization 
do not significantly affect the result, but the width of the a is not small. 

A popular approach to the problem is based on the so-called inverse ampli- 
tude method. Applying it to the xPT representation of the scattering amplitude [36) 
invariably produces a pole in the right ball park. The procedure definitely 
improves the quality of the two-loop approximation of xPT on the real axis, 
because it respects unitarity in the elastic region. The extrapolation into the 
complex plane, however, also contains a number of fake singularities. Like all 
other parametrizations, this approach relies on a model. I do not know of a 
way to estimate the systematic uncertainties generated if QCD is replaced by 
one model or the other - if error estimates are given at all, these necessarily 
rely on guesswork. A thorough discussion of the problems inherent in the IAM 
approach appeared several years ago. [37] 

We have found a method that does not require a parametrization of the data 
at all. It relies on the fact that (a) the ^-matrix has a pole on the second sheet 
if and only if it has a zero on the first sheet, (b) the Roy equations are valid 
not only on the real axis, but in a limited domain of the first sheet, (c) the 
poles from the lowest resonances, a, p, /o(980), all occur in that domain. The 
numerical evaluation of the pole position is straightforward. The one closest to 
the origin occurs at [5] 

M a - % - T a = Vi7 = 441 ^ 8 6 - i 272 ^ 9 2-5 MeV , (9) 

where the error accounts for all sources of uncertainty. We may, for instance, 
replace our representation for S° by the one proposed by Bugg|38j or replace 
the entire scattering amplitude by the parametrization in KPY.J3T] In either 
case, the outcome for the pole position is in the above range. For more details, 
I refer to CCL[2] and to a recent conference report. [39] 

In the meantime, the method described above was applied to the case of irK 
scattering, with the result that the lowest resonance in that channel occurs at 
m K = (658 ± 13) - i (278.5 ± 12) MeV.gD] Evidently, the physics of the k is very 
similar to the one of the a. 



7 On the working bench 

We are currently extending the work described in CGL to higher energies. The 
price to pay is that the contributions from the high energy region then become 
more important. Since the first few terms of the partial wave expansion do not 
represent a decent approximation there, one instead uses a Regge representation 
for the asymptotic domain. In CGL, we borrowed that from Pennington. |41j In 
the meantime, we have performed a new Regge analysis, invoking experimen- 
tal information as well as sum rules to pin down the residue functions. Brief 
accounts of this work were given elsewhere. |42) Figure [3] shows the preliminary 
results obtained for the total cross sections with t-channel isospin and 1: 

^ol = l A + °lot + ¥tot , "III = Wtot + ¥\ot - Hot ■ (10) 
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Figure 3: Asymptotic behaviour of the total cross sections with I t = and 



We include pre-asymptotic contributions and assume that the Regge parametriza- 
tion yields a decent approximation above 1.7 GeV. At lower energies, the un- 
certainties in that representation become larger than those in the sum over the 
partial waves, which is also shown. The comparison with the representation 
used in ACGL confirms the cross section with I t — I, while the one for I t — 
comes out larger by 1 or 2 a. In the solutions of the Roy equations, the effects 
generated by such a shift are barely visible below 800 MeV. We hope to com- 
plete the analysis soon, as well as the application to the electromagnetic form 
factor of the pion, for which an accurate representation is needed in connection 
with the Standard Model prediction for the magnetic moment of the muon. 

It is a pleasure to thank the organizers of the meeting for their kind hospi- 
tality during a very pleasant stay at Chapel Hill and Claude Bernard, Leonardo 
Giusti and Steve Sharpe for useful comments about the text. Also, I acknowl- 
edge Balasubramanian Ananthanarayan, Irinel Caprini, Gilberto Colangelo and 
Jiirg Gasser for a most enjoyable and fruitful collaboration - the present report 
is based on our common work. 
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